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We study the well formulation of the initial value problem of /(_R)-gravity in the metric-affine formalism. 
The problem is discussed in vacuo and in presence of perfect-fluid matter, Klein-Gordon and Yang-Mills fields. 
Adopting Gaussian normal coordinates, it can be shown that the problem is always well-formulated. Our results 
refute some criticisms to the viability of /(i?) -gravity recently appeared in literature. 
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(SI , I- INTRODUCTION 

U 

Every theory of physics is "physically" viable if an appropriate initial value problem is suitably formulated. This means 

■ that, starting from the assignment of suitable initial data, the subsequent dynamical evolution of the physical system is uniquely 
| determined. In this case, the problem is said well-formulated. For example, in classical mechanics, given the initial positions 

■ and velocities of the particles (or of the constituents of a physical system), if the system evolves without external interferences, 
the dynamical evolution is determined. This is true also for field theories as, for example, the Electromagnetism. However, also 
if the initial value problem is well-formulated, we need other further properties that a viable theory has to satisfy. First of all, 
small changes and perturbations in the initial data have to produce small perturbations in the subsequent dynamics over all the 



space-time where it is defined. This means that the theory should be "stable" in order to be "predictive". Besides, changes in 
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the initial data region have to preserve the causal structure of the theory. If both these requirements are satisfied, the initial value 
problem of the theory is also well-posed. 

It can be shown that General Relativity has a well-formulated and well-posed initial value problem but, as for other relativistic 
field theories, we need initial value constraints and gauge choices in order to make Einstein's field equations suitable to correctly 
formulate the Cauchy problem. The consequence of this well-posedness is that General Relativity is a "stable" theory with a 
robust causal structure where singularities can be classified. A detailed discussion of these topics can be found in IHIH]- 
, In this paper, we focus attention on the well formulation of the initial value problem of metric-affine /(i?)-theories of gravity, 
("T) 1 where f(R) is a generic function of the Ricci scalar R. The aim is to prove that the Cauchy problem is well formulated in 
* \ different and physically important cases. 

Thepossible modifications of Einstein's theory has a long history which reaches back to the early times of General Relativity 
fiSlaSfe, 9]. The early extensions and modifications were aimed at a unification of gravity with other fundamental interactions 
. (in particular the Electromagnetism). The recent interest in such modifications comes from cosmological observations. For 
comprehensive reviews, see for example ifToL [Til [l2l [l3ll . These observations usually lead to the introduction of additional ad- 
hoc concepts like dark energy/matter, if interpreted within Einstein's theory. On the other hand, the emergence of such stopgap 
measures in a cosmological context could be interpreted as a first signal of a breakdown of General Relativity on these scales 
OH EH, and led to the proposal of many alternative modifications of the underlying gravity theory. In particular, more recent 
works focused on the cosmological implications of / ( .R W ravit y since such models may lead to an alternative explanation of 
the acceleration effect observed in cosmology lfl6l n% [T8l Il9i 12011 and to the explanation of the missing matter puzzle observed 
at astrophysical scales |fl5l , l2lll . 

While it is very natural to extend Einstein's gravity to theories with additional geometric degrees of freedom, see for example 
Il22ll23ll24ll for some general surveys on this subject as well as 1125m 1 for a list of works in a cosmological context, recent attempts 
focused on the old idea of modifying the gravitational Lagrangian in a purely metric framework, leading to higher-order field 
equations. Due to the increased complexity of the field equations in this framework, the main body of works dealt with some 
formally equivalent theories, in which a reduction of the order of the field equations was achieved by considering the metric and 
the connection as independent objects IB^l . In addition, many authors exploited the formal relationship to scalar-tensor theories 
to make some statements about the weak field regime, which was already worked out for scalar-tensor theories 112711 . 

However, a concern which comes with generic higher-order gravity theory is linked to the initial value problem. It is, up 
to now, unclear if the prolongation of standard methods can be used in order to tackle this problem in every theory. Hence it 
is doubtful that the full Cauchy problem can be properly addressed, if one takes into account the results already obtained in 
fourth-order theories stemming from a quadratic Lagrangian j28l l29ll . 

On the other hand, being /(i?) -gravity, like General Relativity, a gauge theory, the initial value formulation depends on 
suitable constraints and on suitable "gauge choices" that mean a choice of coordinates so that the Cauchy problem results well- 
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formulated and possibly well-posed. In [28, 30], the initial value problem was studied for quadratic Lagrangians in the metric 
approach with the conclusion that it is well-posed. On the other hand, in 03111 . the Cauchy problem for generic /(i?) -models 
has been studied in metric and Palatini approaches using the dynamical equivalence between these theories and the Brans-Dicke 
gravity. The result was that the problem is well-formulated for metric approach in presence of matter and well-posed in vacuo. 
For the Palatini approach, instead, the Chauchy problem is not well-formulated even in vacuo since, considering the 3 + 1 ADM 
formulation of the equivalent scalar-tensor theory, the equivalent Brans-Dicke parameter oj = —3/2 leads to a non-dynamical 
equivalent field <f> and to the impossibility of a first-order formulation of the problem. 

In this paper, we will adopt a different approach by which it is possible to show that the Cauchy problem of metric-affine 
/(i?)-gravity is well-formulated and well-posed in vacuo, while it can be, at least, well-formulated for various form of matter 
fields like perfect-fluids, Klein-Gordon and Yang-Mills fields. Therefore, our results refute the conclusions stated in 13111 on the 
viability of metric-affine /(i?)-theories of gravity. The reason of the apparent contradiction with respect to the results in ll3lll 
lies on the above mentioned gauge choice. Following |2J], we adopt Gaussian normal coordinates. Such a choice, introducing 
further constraints on the Cauchy data surface, results more suitable to set the initial value problem in such a way that the well- 
formulation is always possible. Of course, in order to prove the complete viability of metric-affine /(i?) -theories of gravity, 
the well posedness of the Cauchy problem has to be shown too. This topic is not dealt with in the present paper, being still an 
open problem under investigation. Anyway, here the important point is that, up to now and in our opinion, no objections to the 
viability of metric-affine /(i?) -gravity based on the ill-formulation of the initial value problem can be made. 

The layout of the paper is the following. In Sec. II, following J2l, the initial value formulation of General Relativity is recalled. 
Sec. Ill is devoted to set the problem for /(i?)-gravity in metric-affine formalism in vacuo. In this case, the problem is well- 
formulated and well-posed since dynamics reduces to the Einstein theory plus cosmological constant, /(i?) -gravity in presence 
of matter is discussed in Sec. IV. As paradigmatic examples we discuss the coupling with perfect-fluid matter and Yang-Mills 
fields, showing that in both the cases the initial value problem results well-formulated. 

In Sec.V, the problem in presence of a Klein-Gordon scalar field is discussed. In all these cases, the Gaussian normal co- 
ordinates allow the well-formulation of the initial value problem. It is important to stress that further constraint equations can 
emerge on the initial surface and this fact could reduce the set of admissible initial data. Discussion and conclusions are given in 
Sec.VI. The Appendix A is devoted to the demonstration of two useful propositions related to the effective stress-energy tensor 
and the Bianchi Identities. 



II. WELL FORMULATION OF THE CAUCHY PROBLEM FOR GENERAL RELATIVITY 

Before starting with our considerations for /(i?)-gravity, let us recall the initial value formulation of General Relativity (i.e. 
f(R) = R) where it is well-formulated (and also well-posed as shown in [3]). We adopt the formalism developed in |2|. 

Let us consider a system of Gaussian normal coordinates (3[] where the Latin indexes i, j run from 1 to 4 and the Greek indexes 
a run from 1 to 3. In these coordinates, the time components of metric tensor are 1744 = — 1 and g 4a = with the signature 
(+ + H — ). These are particularly useful to split the spatial hypersurface S3 from the orthogonal time-geodesics in a given 
space-time M, 

After, given a second rank symmetric tensor Wij, defined on the globally hyperbolic space-time manifold (M,gij), it is 
possible to define the symmetric conjugate tensor W*j as W*j = Wij — \Wgij, where W := W^ g^ is the trace of Wij. 
Furthermore, if Si is a space-time domain in M where 344 ^ and S3 is the 3-surface given by the equation x 4 = 0, then the 
following statements are equivalent: 



a) 


Wij: 


= in Si. 




b) 




= and W 4j 


= in Si. 







= and W] v 


= in S 4 with W 4j = in S 3 



Let us take into account now the Einstein equations in the form 

dj = -fcly with the Bianchi indentities T*f = , (1) 

where Gij — Rij — —Q%jR is the Einstein tensor. We can define the tensor 

Wij : Gij + kTij . (2) 

The conjugate tensor is 



W* = Rij + kT* , (3) 
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and then the Einstein equations becomes 

Wij = . (4) 

These are 10 equations for 20 unknown functions g^ and . Let us assign now the 10 functions ^4 and T a p. The remaining 10 
functions g a p and Ta, can be determined by Eqs. (|4j. Using the above results, these functions can be expressed in the equivalent 
form 

R a p + kT* = O, W^=T; {i = 0, (5) 

with the condition 

+ kT^j = on the surface x 4 — . (6) 

Eqs. © can be rewritten in the form 

5a/3,44 = 2R a/3 - ^Ag a/3A + g^ u g a ^ A 90uA + 2kT* afi , (7a) 

I ' ; ' .... • 1 ... 1 j _ 1 lj J a; ( /b ) 

where _R Qi 3 is the intrinsic Ricci tensor defined on the surface x 4 — 0, T i4 k is the Levi-Civita connection related to the metric 
gij, the coma is the partial derivative and 

A := g^g^i . (8) 

In the same way, the constraint equations © become 

A, a - D a g aaA + 2kT 4a - , (9a) 

R- ^A 2 + jB + 2kT A4 = 0, (9b) 

where R is the intrinsic curvature scalar of the surface x 4 = 0, D a denotes the covariant derivatives on the surface x A = 
associated to the Levi-Civita connection of the intrinsic metric g a (3\ x 4 =0 and 

B = g» v g pa g W Ag V oA- (10) 
Let us assign now the set of Cauchy data on the surface x 4 = 

ga/3, gaf3A, T i4 . (11) 

Such data have to satisfy the constraint equations Eqs. (0 explicitly give the values of the quantities 

5ct/3,44, ^4j,4 , (12) 

as a function of the Cauchy data. By deriving Eqs. (0, it is straightforward to obtain the time-derivatives of further order as a 
function of the Cauchy data. This procedure allows to reconstruct the solution of the field equations as a power-law series of the 
time variable x A . 

In other words, this means that the 3-surface S3, given by the equation x A = 0, is a Cauchy surface of the globally hyperbolic 
space-time (M,gij) and that the initial value formulation (the Cauchy problem) is well-formulated in General Relativity. Our 
task is now to extend these results to /(i?)-gravity in metric-affine formalism. 



III. THE CAUCHY PROBLEM FOR f(7?)-GRAVITY IN METRIC-AFFINE FORMALISM IN EMPTY SPACE 

In the metric-affine formulation of /(i?) -gravity, the dynamical fields are given by the couple of functions (g, V) where g is 
the metric and V is the linear connection. In vacuo, the field equations are obtained by varying with respect to the metric and the 
connection the following action 



A(g,T) = / y/\g\f(R)ds (13) 
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where f(R) is a real function, R(g,T) = g 1 ^ Rij (with Rij := R h ih j) is the scalar curvature associated to the dynamical 
connection T. 

More precisely, in the approach with torsion, one can ask for a metric connection T with torsion different from zero while, in 
the Palatini approach, the T is non-metric but torsion is null B32I1 . 

In vacuo, the field equations for /(i?) -gravity with torsion are lf32L [33l l34ll 

f\R)Ri3-lf(R)aa = o, (14a) 



2/' dxP 

while the field equations for /(i?)-gravity a la Palatini are I26ll35l,l36l l3 



T h = (S p S h - 5 p S h ) (14b) 



f{R)Ra - ^f{R)g l3 = o , (15a) 



V k (f'(R)g ij ) = 0. (15b) 

In both cases, considering the trace of Einstein-like field equations ( 114at e (|15a| i, one gets 

f'(R)R-2f(R) = 0. (16) 

It is easy to conclude that scalar curvature R is a constant coinciding with the solution of Eq.([T6ll. In this case, Eqs. ( 1 14-bb e (I15bl ) 
imply that both dynamical connections coincide with the Levi-Civita connection associated to the metric g,j which is solution of 
the field equations. In other words, both theories reduce to the Einstein theory plus cosmological constant. As it is well known, 
in this case the Cauchy problem is well formulated and well posed too |2|]. 



IV. THE CAUCHY PROBLEM IN PRESENCE OF MATTER 

Let us now take into account the presence of perfect-fluid matter in the formulation of the Cauchy problem for f(R) -gravity. 
In order to deal simultaneously with the Palatini approach and the torsion, we will consider the connection not coupled with 
matter. In other words, we will assume that the matter Lagrangian does not explicitly depend on the dynamical connection. With 
this working hypothesis, the field equations are 

f'(R)Ri3 - \f{R)9ij = , (17a) 



with 



in the case of /(i?)-gravity with torsion, and 



f'(R)Rij - nW)9ij = Ey , (18a) 



V k (f(R) 9ij ) = 0, (18b) 

1 SjO 

in the case of /(i?)-gravity in the Palatini approach. In Eqs. (117at and (|18a| i, Sj» := — m is the stress-energy tensor. 

Considering the trace of Eqs. (I17ab and (1 1 8 at . we obtain a relation between the curvature scalar R and the trace of the stress- 
energy tensor S := g^ Sy. We have 

f'(R)R-2f(R) = X. (19) 

It is worth noticing that any time that E = const, the theory reduces to the General Relativity with cosmological constant and 
the initial value problem is identical to the above empty-space case. 
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By the hypotheses that the relation $1% is invertible and E ^ const., the curvature scalar R can be expressed as a function of 
E, that is 

R = F(E) . (20) 

Starting from this fact, it is easy to show that the Einstein-like equations of both Palatini and metric-affine theory with torsion 
can be expressed in the same form I32ll33ll37ll . that is 



1b„ _ 1 ^ , 1 ( 3 dl P d( P , dl P , 3 9( P dl f hi- 



(21) 



where we have defined the effective potential 



1 



V(<p) := - [pF-^if')- 1 ^)) + V 2 (f'rH<P)] , (22) 
for the scalar field 

<p := /'(F(E)) . (23) 
Introducing the conformal transformation g^ = (pgij, Eq. d2TT > assume the simpler form (see for example ll32i[37l[38ll ) 

Rij - \R9H = ~£y - ^Vi^Sij > (24) 

where R\j and R are respectively the Ricci tensor and the curvature scalar derived from the conformal metric gij . It is worth 
noticing that the conformal transformation is working if the trace E of the stress-energy tensor is independent of the metric g^ . 
In such a case, Eqs.d24l> depend exclusively on the conformal metric g^ and the other matter fields. 

Furthermore, we have to stress that the connection V, solution of the field equations with torsion, is given by 



r h -t h + —^S h -—^aP h a- (25) 



where is the Levi-Civita connection induced by the metric g^, while the connection T, solution of the Palatini field equa- 
tions, coincides with the Levi-Civita connection associated to the conformal metric </y . Both connections, T and T, satisfy the 



relation 



The identities 



^ = T » + 2^^- (26) 



^ 1 « + 2^ dxi 1 2^ dxP 9 9lJ + 2 V dx> { ' 

hold. This result shows the relation between the Levi-Civita connections induced by the metrics g^ and g^ . 

In general, the Einstein-like Eqs. (fSTJ have to be considered together with the matter field equations. To this purpose, we have 
to keep in mind that the conservation equations for both the metric-affine theories (with torsion and a la Palatini) coincide with 
the standard conservation laws of General Relativity (see the Appendix A). This means 

Vj>: ij = 0. (28) 

It is straightforward to show that Eqs. d28l) are equivalent to the conservation laws 

y 3 T^=0 where T« = -E 4J ■- -^V(cp)g i3 , (29) 

for the conformally transformed theories (l24t . In fact, by an explicit calculation of the divergence WjT 1 ^ where the relations 
d27| i have been used, we obtain the equations 

VT V = ^E 4 , + (-iE + iv(„) - Vfc)) . (30) 
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The constraint equations (f28b and d29i l are then mathematically equivalent in view of the relation 

V--V(<p)+2V'(<p) =Q, (31) 

V 

which is equivalent to the definition ip = f'(F(Jl)) ll32ll . 

With these results in mind, the Cauchy problem for Eqs. (|2TT > and the related matter field equations can be faced by discussing 
the equivalent initial value problem of the conformally transformed theories. Considering as in the above General Relativity 
case, Gaussian normal coordinates and starting from Eqs. J24b and ( [29) , it is easy to conclude that the Cauchy problem is well 
formulated also in this case. 

It is worth pointing out that, in general, the matter field equations imply the Levi-Civita connection derived from the metric 
gij and not the connection induced from the conformal metric g^. In view of the relation (f27), this is not a problem since the 
connection T can be expressed as a function of the connection f and of the scalar field ip, which, on the other hand, is a function 
of the source matter fields. As a result, we could obtain slightly more complicated equations implying further constraints on the 
initial data but, in any case, the same equations can be always recast in "normal form" with respect to the maximal order time 
derivatives of matter fields allowing a well-formulated Cauchy problem fU. 

As an example, let us examine in detail the perfect-fluid case with equation of state p = p(p). The corresponding stress-energy 
tensor is 



E« = PViVj + P (VM + gii ) , (32) 
the matter field equations are given by Eqs. d28l l with the further condition 

.</,.. I- (33) 

where V 3 are 4-velocities. Specifically, Eqs. ( |28l l give the field equations 



(p + pV 3 ) ! . Vi + (p + p) V i{j V 3 + ^=0, (34) 

where | denotes the covariant derivative with respect to the Levi-Civita connection induced by g^. By saturating with V 1 , one 
gets 

(pV 3 ) { . = -pV 3 j , (35a) 
while by substituting Eq. (I35ab into Eqs. ( 134b for a = 1, 2, 3, we have 

(p + p) V j V{* = -|| (V a V 3 + g<* 3 ) . (35b) 

As already pointed out, Eqs. d33l and d35b involve the metric g^ and its first derivatives. By using the relations d27b . we can 
rewrite these equations in the conformal metric g^, the scalar field function (of the matter density) = (f(p) and their first 
derivatives. Immediately we get 

- 9lJ V l Vi = -l, (36a) 
f 



(py j ) 



2 dip 
ip dx s 



pv s 



(d\P_ 



T^V 3 



2 dip 
ip dx s 



V 



(36b) 



(p + p) V 



dV a 
dxi 



- jf 



'V s 



1 

2~p 



9f 5 a 

dx s 3 



dip 



5? 



dip 
dxP 



pot 

9 9js 



V s 



.I*L (v a v j 



(36c) 



±1 (depending on the sign of ip) and g a 4 — 0, we obtain, from Eq. (|36a| ), 
Eqs. ( I36ab and (I36bb can be considered as linear equations for the 



In the Gaussian normal coordinates, where 1744 
the expression of V 4 in terms of the remaining V a 

dV a dp 

functions . and T — T , The explicit resolution of these equations, in terms of the unknown functions, could give rise to 

dx^ OT 4 

further constraints on the initial data and on the form of the /(i?) -function. In any case, from Eqs. ( I36ab and ( I36bb and in 
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Gaussian normal coordinates, one can derive the quantities — 4 e — j as functions of the initial data g a p, g£f , and p 

allowing to put the matter-field equations in normal form. This means that the Cauchy problem is well-formulated. 

As another important example, let us consider the initial value formulation of /(i?)-gravity coupled with Yang-Mills fields, 
in particular with Electromagnetic field. Also in this case, the problem is well-formulated. In fact, the stress-energy tensor of 
a Yang-Mills field has null trace. From Eq.(|T9b, it is easy to prove that the curvature scalar R is constant and then, by using 
Eqs. ( 117bt and ( 118b| ), it is straightforward to conclude that the dynamical connection coincides with the Levi-Civita connection 
of the metric </y . In this situation, both theories (a la Palatini and with torsion) reduce to the Einstein theory with cosmological 
constant and then the Cauchy problem is well-formulated (this last result is immediate for the Electromagnetic field). Regarding 
the well-posed initial formulation, the above results work for any theory where the trace of stress-energy tensor is a constant 
since, as above, the /(i?)-gravity reduces to Einstein gravity plus cosmological constant. 



V. THE CAUCHY PROBLEM IN THE CASE OF COUPLING WITH A SCALAR FIELD 

Let us consider now the case of coupling with a Klein-Gordon scalar field ip with self-interacting potential U (tp) — —m 2 ip 2 . 
The stress-energy tensor is given by 

dip dip 1 ■ { dtp dip 2 2 \ 

= Tt — - "7: — r q J — — — — q l H + m ip . (37) 

The trace of tensor ( f37T > is 

dxP dxi 

In the relation (138b . the metric g pq explicitly appears. This means that the conformal transformation procedure of the previous 
section cannot be applied in this case and then the well-formulation of the Cauchy problem has to be directly shown starting 
from the field Eqs. (L?Th and the Klein-Gordon equation (see also [40]), that is 



As standard, we take into account Gaussian normal coordinates where g± a = and 344 = — 1. It is easy to show that = 
and then Eq. d39l ), suitably developed, can be rewritten in the form 

- .9 ' [ a „„ fl „ fl ~ r a/3 TOT - m V- , (40) 



(5a: 4 ) 2 \dx a dxP 

where 



- h chp_ = 1 Xy ( dg ai dg Pl _ dg af3 \ dip 1 dg af3 dip 
af3 dx h 2 9 \ dx? dx a dx"/ J dx x 2 5a; 4 dx 4 ' ' 

The Einstein-like equations are exactly of the form (fSTJ but now they depend also on tp 

1 6„ _ X ^ , 1 ( Z dtp dtp ^ ^ dtp t Z dtp dtp ^ hk 

lij 

(42) 



Rij ~ 2 R 9v ~ -S« + — 2 \-2dx^dxl + ^dtf + Idx^dx^ 9 9ij 



'^dx* 913 ~ V ^ 9ij 



Adopting the notation of Sec. II, we can rewrite Eqs. d42b in the form 

Wij := tpG i:l + Tij = , (43) 

where 



_ 1 / 3 dtp dtp ~ dp 3 dtp dtp hk 

13 ~ 13 + Tp \2dx~i~dxl + ^'dx 1 + Zdl^dx* 9 913 



(44) 
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We have now all the ingredients to discuss the Cauchy problem for the system d39l l and (143V As first result, it has to be pointed 
out that Eqs. fl39l > imply the conservation laws 

Vj£ i3 '=0. (45) 
In Appendix A, it is shown that Eqs. (PES! are equivalent to the conservation laws 

V j Wij = 0, (46) 



see also 13911 . In this case, it is useful to remember that both metric-affine theories (a la Palatini and with torsion) are equivalent 
to a metric scalar-tensor theory with a Brans-Dicke parameter u = — § IT321 l35l l36Tl . 

For the the Cauchy problem, it is important to point out that the second of Eqs. (0 is given by the Klein-Gordon Eq. (f39b - 
Regarding the constraint ©, from the field Eqs. (142b . it is easy to show that it involves only first order derivatives with respect to 
the time variable x 4 . In order to show this point, let us observe that the r.h.s. of Eqs. d42l) contains the second derivatives of the 
scalar field ip = f'(F(E)) and then implies the second derivatives of the metric gij and the third derivatives of the Klein-Gordon 
field tp. In more details, the second derivatives of ip are present in the terms 

- gy = ay _ f s 3y_ 

3 8x l dx^dx 1 31 3x s ' 

and 

V ^=^^-^ r -^)^" (4?b) 

Assuming as before Gaussian normal coordinates, which means 544 = g 44 = — 1 and g 4ce = g 4a = 0, it is easy to show 

d 2 ip dip 
that in the equations Ga = — T 4 .; there are no second order time derivatives . , and the only time derivatives are r 

(ox^y ~ ox 4 

and — — 7— — . On the other hand, it is useful to observe that these terms involve also second order time derivatives of the 

dx 4 dx a 

Klein-Gordon field ,„ ,^„ , and — — — ,.„ , the second order spatial derivatives of the time derivative — — - „ - . , and the 

{dx 4 ) 2 dx a {dx 4 ) 2 F dx a dxf>dx 4 

a ■ ft u >■ d 9 a P a d2 9aP 

derivatives of the metric - „ and 



hh d( P „ ( 'hk 9 V hkr s 9(p 



dx 4 dx-rdx 4 ' 

At this point, it is sufficient to observe that the Klein-Gordon Eq. d40l i, evaluated on the initial surface x 4 — 0, allows to 

d 2 ib d 3 ip 
express the second order time derivative — - , NO and its spatial derivatives — — j-r on the same surface as a function of the 

F {dx 4 ) 2 F dx a {dx 4 ) 2 

Cauchy data (g a/3 , ip, J^j) and of their spatial derivatives. The remaining quantities g a gjgg 4 » and dx~>dx 4 ' 

defined on x 4 = 0, can be directly calculated as functions of the same Cauchy data. 

In conclusion, the equations Ga — — on x 4 = involve only the Cauchy data and their spatial derivatives. In other words, 
they can be considered as genuine constraint equations for the Cauchy data. 

At this point, some further remarks are necessary regarding the corresponding Eqs. (O for the case considered. We have 
already seen that Eqs. d42b imply second derivatives of the metric and third derivatives of the scalar field if). Our task is 
now to prove that Eqs. © imply at most second order time derivatives (i.e. with respect to the variable a; 4 ) of the metric and at 
most first order time derivatives of the Klein-Gordon field. The first request is obvious since only second order derivatives of the 
metric are involved. Besides, the second request is satisfied thanks to the Klein-Gordon Eq. ( l40b . In fact, by deriving Eq. (f40b 
with respect to the variable x 4 and with respect to the remaining variables x a , by substituting Eq. (l40t into the obtained results 

3 ib 3 ib 

and evaluating all the quantities on the initial surface x 4 = 0, it is possible to express — and — — — on x 4 = as 

(dx 4 ) 3 dx a {dx 4 ) 2 

d 2 g a 

functions of the Cauchy data and, at most, of the derivatives .[ . 

d{x i Y 

In conclusion, in the equations corresponding to Eqs. (0, only second order time derivatives of the metric appear. Obviously, 

d 2 g 

the initial value problem is well-formulated if these equations are in normal form with respect to at least on x 4 = 0. 

<9(x 4 ) 2 

As in the perfect-fluid case, such a request could impose further constraints on the initial data and, possibly, on the form of the 
/(i?)-function, but, in general, it works and the corresponding Cauchy problem is well-formulated. 
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VI. DISCUSSION AND CONCLUSIONS 



In this paper, we have shown that the initial value problem for metric-affine /(i?)-gravity is well-formulated. This means 
that there are no objections to the viability of metric-affine /(i?) -theories of gravity based on the well formulation of the initial 
value problem. However, also the well-posed problem is necessary in order to achieve a complete control of dynamics but, in 
this case, the role of source fields has to be carefully discussed. This topic will be discussed in a furthcoming paper. 

Since /(i?)-gravity is a gauge theory, like General Relativity, it is crucial the choice of suitable coordinates in order to correctly 
formulate the problem. We have adopted Gaussian normal coordinates which can be defined any time that the derivative operator 
Vj arises from a metric </y. They are also called synchronous coordinates and are particularly useful for calculations on a given 
non-null surface S3, i.e. an 3-dimensional embedded sub-manifold of the 4-dimensional manifold M. They allow to define 
uniquely orthogonal geodesies to S3 and then to correctly formulate the conditions of validity for the Cauchy-Kowalewski 
theorem . 

In metric-affine formalism, it is always possible to show that a given /(i?)-theory, in vacuo, can be recast in Einstein's gravity 
plus a cosmological constant. This means that the initial value problem is always well-formulated and well-posed. The same 
conclusion holds in the case of matter coupling every time that the trace of the stress-energy tensor is a constant. 

As shown in ll32[[33ll34ll . by introducing matter fields in the Palatini and in the metric-affine approach with torsion, one can 
define, in two steps R = F(E) and p :— /'(F(E), a suitable scalar field that allows: i) to reduce the theory to scalar-tensor 
theory; ii) to relate the form of f(R) to the trace of the matter (field) stress-energy tensor. In this case, it is always possible a 
well-formulation of the Cauchy problem avoiding the singularities which could emerge with other gauge choices 13111 . Besides, 
matter fields could induce further constraints on the Cauchy surface x A which, if suitably defined, lead to the normal form 
of the matter-field equations. This is one of the main requests to obtain a well-formulated value problem. However different 
fields, acting as sources of the field equations, like perfect fluids, Yang-Mills, and Klein-Gordon fields, could generate different 
constraints on x A . Such constraints could imply also restrictions on the possible forms of f(R). In conclusion, as in General 
Relativity, the gauge choice is essential for a correct formulation of the initial value problem while the source fields have to be 
carefully discussed to obtain a well-posed problem. In a forthcoming paper, we will study specific problems and /(i?) -models 
where both the issues (i.e. well-formulation and well-posedness) could be achieved. 

APPENDIX A: THE EFFECTIVE STRESS-ENERGY TENSOR AND THE CONSERVATION LAWS 



Proposition A. 1. Eqs. ( 12 U . ( 122b and (1231 > imply the usual conservation laws V-'Sy = 
Proof. First of all, we recall that Eq. (l22l and (l23l is equivalent to the relation 

£ - -VUp) + 2V'(<p) = 

(see i32ll for the proof). After that, taking the trace of Eq. (fJTJ into account, we get 



(Al) 



E = -<pR - l-<Pi<P* + 3V<^ + -V(ip) 

Zip. if 

Q 

where for simplicity we have defined pi :— ■ 7 —. Substituting Eq. (IA2t in Eq. (IAU . we obtain 

ox 1 



(A2) 



R + Wfif* - -v^ + \v(ip) - -v'(ip) = 



2p 2 



(A3) 



We rewrite Eq. (|2TT > in the form 



f 



pRij ~ n R 9i] = s y + - —zfiVj + <pVj<Pi + ~;Vh<P 9ij- 



(A4) 



-pV h p h g i:j - V(p)g lj 



The covariant divergence of dA4b yields 



1 - - 



-V 3 



3 3 

--p t Pj + -p h p h g i:j - V(<p)gij 



(A5) 
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By definition, the Einstein and the Ricci tensors satisfy V-'Gy = and (V J (/?)i?y = (v^'VjVi - ViV^V,-") if. ThenEq. <|A5]> 



reduces to 



Finally, making use of Eq. ( IA3b it is easily seen that 

- -RVnp = V 3 



1 



1 /" 3 3 h T ,, , "\ 

- I —jViVj + ^<PhV 9ij ~ VW)9ij\ 



3 3 

■■^PWj + ~£PW 9%j ~ V(ip)g it 



(A6) 



(A7) 



from which the conclusion V J = follows □ 



Proposition A.2. Given the Levi-Civita connection 

f h -f h + ±^S h - —^aP h a - + —^S h 
associated with the conformal metric tensor g = (pg, and given the effective energy-impulse tensor 

Tij = -E« - -^Vi^gij 



(A8) 



(A9) 



the condition V 3 Tij = is equivalent to the condition V^Ey = 



Proof. 



We have separately 



1 



V'7„ = -,r'V,//,.. = -'/-•• 



1 



VJ]l ~ 2~ V [JJ? 6 4 — 



6" 

dx s 1 



^a u «a ■ \ T ■+ 
q x u9 S 'si M<?j + 



1 



-//'~'v.r, ( = -V'v.. ( -v,.,. _ v { ,- )ttl! j = _ v-'s« - 



i 



i . 



i 



2^ dx s 3 dx 



1 ^ 



•■P 



ip 3 dx 



(A10) 



(All) 



y 9 2tp {dx* s + dx* di dx- 9 9sl ' lqj 



tp 2tp ydx 1 



dx s " 1 
1 



-I s .' 



2^ 



9a; 

- -V(<p)g S j ^ 



dx 



dx t 



dx s 9.x 



1 ^ s 2 

2^ 3 dx* dx 1 W) 



(A12) 



I s .' 



tp 2ip \dx 3 



tp 2ip \dxi 



dx s 3 
1 
2~ 



dx 1 



9 uq 9s 3 



dx s 3 
1 



S i9 2^(^)5* 

V If 



dtp 



dtp 



(A13) 
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Collecting Eqs. (IA11I ), ( IA12b and ( IA13I ) we have then 



2 95 



(A14) 



because the identity — i£ + — V(<y3) — V^'(ip) = holds identically, being equivalent to the definition ip = f'(F(T,)) [32]. □ 
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